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Abstract 

We present a multi-lattiee kinetie Monte Carlo 
(kMC) approaeh that effieiently deseribes the 
atomistic dynamics of morphological transitions 
between commensurate structures at crystal sur¬ 
faces. As an example we study the reduction of 
a (\/5 X y/5)R21° PdO(lOl) overlayer on Pd(lOO) 
in a CO atmosphere. Extensive density-functional 
theory calculations are used to establish an atom¬ 
istic pathway for the oxide reduction process. 
First-principles multi-lattice kMC simulations on 
the basis of this pathway fully reproduce the ex¬ 
perimental temperature dependence of the reduc¬ 
tion rate [Fernandes eta/., Surf. Sci. 2014 , <527, 
31-39] and highlight the crucial role of elementary 
processes special to the boundary between oxide 
and metal domains. 

1 Introduction 

In recent years first-principles (Ip) kinetic Monte 
Carlo (kMC) simulations have established a new 
standard for the microkinetic modeling of sur¬ 
face chemical processes and heterogeneous catal¬ 
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ysis.-v^ On the one hand, this approach employs 
material-specific rate constants for the elemen¬ 
tary processes, presently typically computed by 
density-functional theory (DFT) and transition- 
state theory (TST). On the other hand, and in 
contrast to still prevalent mean-field rate equation 
based microkinetic models, the Ip-kMC approach 
explicitly accounts for spatial heterogeneities, cor¬ 
relations, and fluctuations of the atomic structure. 
The downside of the approach is the computational 
cost of the underlying DFT-TST calculations. In 
practice, Ip-kMC simulations are therefore hith¬ 
erto performed on lattice models. This exploits the 
translational symmetry of the crystalline lattice to 
dramatically reduce the required number of differ¬ 
ent Ip rate constants.-^ For simple reaction net¬ 
works like the CO oxidation reaction at low-index 
single-crystal model catalyst surfaces, this num¬ 
ber of Ip rate constants can be as low as 10-20, 
comprising adsorption, desorption, diffusion and 
reaction steps at and between the different lattice 
sites.— 

By construction, the conventional lattice Ip- 
kMC approach cannot directly be applied to (sur¬ 
face) morphological rearrangements, i.e. any type 
of phase transition that includes changes in the 
considered lattice site arrangement. Notwithstand¬ 
ing, for the (late) transition-metal catalysts typi¬ 
cally employed in oxidation catalysis, a possible 
formation of (surface) oxides has been discussed. 
Possibly, especially the continuous formation and 
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reduction of such oxidic films could be a major ac¬ 
tuator for the observed catalytic activity.—-— As 
a first step to make such systems and phenomena 
accessible to Ip-kMC we here describe a multi¬ 
lattice Ip-kMC approach, which is applicable to 
morphological transitions between commensurate 
crystalline lattices. The simple concept is that we 
work on a super-lattice that contains both lattices 
simultaneously. In a spatial region describing one 
stable phase, only the sites of the corresponding 
lattice are "enabled". The actual morphological 
transition is then described through the elementary 
processes of its atomistic pathway, which "dis¬ 
able" sites of one lattice and "enable" sites of the 
other lattice. 

We illustrate this multi-lattice Ip-kMC idea by 
modeling the surface oxide reduction at Pd(lOO). 
Recent experimental work— observed significant 
changes in the rate of oxide reduction in the tem¬ 
perature range from 300-400 K, when exposing an 
initially prepared (\/5 x y/5)R21° PdO(lOl) sur¬ 
face oxide film (henceforth coined \/5-oxide for 
brevity)—!^ to a 5 x 10“^^ bar CO atmosphere. 
Analysis of the data using fitted mean-field kinetic 
models suggested that the reduction is controlled 
by boundary processes between surface oxide and 
metal domains.— In the present paper, we use ex¬ 
tensive DFT calculations to establish an atomistic 
pathway for the oxide reduction process. Based 
on this pathway, our multi-lattice Ip-kMC simu¬ 
lations perfectly reproduce the observed trends in 
the reduction rates, and show that the critical steps 
in the reduction of the surface oxide correspond in¬ 
deed to cross-reactions between the metal domains 
and the surface oxide. 

2 Theory 

2.1 Kinetic Monte Carlo simulations 
for heterogeneous catalysis 

We focus on thermally driven surface chemical re¬ 
actions that involve a site-specific (covalent-type) 
binding of the reaction intermediates. Such reac¬ 
tions generally follow a so-called rare-event time 
evolution. That is, on time scales that can be simu¬ 
lated by ab initio molecular dynamics all reaction 
intermediates only vibrate around the lattice sites 


to which they are bound. Only rarely, i.e. after 
very many such vibrations, an elementary process 
changes this population at the lattice sites, where 
the duration of this process itself is very short. 
A diffusion process brings a reaction intermediate 
to another site, adsorption and desorption fill and 
empty sites at the surface, surface reactions con¬ 
vert reaction intermediates into another etc. Dur¬ 
ing the long time intervals between such individ¬ 
ual elementary processes the reaction intermedi¬ 
ates equilibrate fully and thereby are assumed to 
forget through which elementary process they ac¬ 
tually reached their current site or state. 

Viewing the entire population at all sites as one 
possible configuration, individual events in form 
of elementary processes thus advance the system’s 
state from one configuration to another. In each 
configuration reaction intermediates are assumed 
to have forgotten their past. In this Markov ap¬ 
proximation events will occur uncorrelated from 
previous ones and the time evolution of the system 
is described by a Markovian master equation 

Pu{t)=Y^WuvPv{t)-Wyupu{t) , ( 1 ) 

V 

where puii) is the probability for the system to be 
in configuration u at time t, and w^u is the transi¬ 
tion rate (in units of time at which configura¬ 
tion u changes to configuration v. For heteroge¬ 
neous catalysis, solving this equation is in prac¬ 
tice not feasible by any direct approach.-*^ The 
kMC approach addresses this problem by gener¬ 
ating configuration-to-configuration trajectories in 
such a way that an average over an ensemble of 
such trajectories yields the correct probability den¬ 
sity pu{t) of eq. [I]'^ Instead of evaluating and 
storing the entire transition matrix w, a kMC code 
thus only focuses on-the-fiy on those transition 
matrix elements Wvu that are actually required to 
propagate the trajectory. 

2.2 Efficient Ip-kMC on one or more 
lattices 

Even with this focus on the required Wyu, a kMC 
simulation can still be an extremely demanding 
task, if the rate constants determining the Wyu are 
to be provided through Ip-calculations. Emerg- 
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ing so-called adaptive or self-learning kMC ap¬ 
proaches—-— address this problem by generat¬ 
ing look-up tables of already computed Wyu in 
the course of the simulation. However, an orders 
of magnitude higher efficiency is needed to make 
comprehensive Ip-kMC simulations tractable for 
systems as complex as the present one. This effi¬ 
ciency is achieved when the sites involved in the 
surface catalysis can be mapped onto a static lat¬ 
tice. Due to the translational symmetry of the lat¬ 
tice, each equivalent site type then features the 
same group of in principle possible elementary 
processes, and these elementary processes typi¬ 
cally depend only on the occupation of nearby 
sites. While this does not reduce the number of 
Wyu required during a Ip-kMC trajectory, it dra¬ 
matically reduces the total number of inequivalent 
Wyu and therewith Ip rate constants that need to 
be computed with DFT-TST. In lattice Ip-kMC 
these Ip rate constants are then in practice precom¬ 
puted and stored, to be looked up on demand in the 
evolving kMC simulation. 

Assuming a static arrangement of sites, lattice 
Ip-kMC can not directly be applied to systems un¬ 
dergoing a morphological transition. In the present 
context, this would be the reduction of the initially 
oxidized surface to a metallic state, where the ox¬ 
ide overlayer and metal substrate exhibit different 
crystal lattices. With our suggested generaliza¬ 
tion to multi-lattice Ip-kMC transitions between 
commensurable lattices become accessible. The 
essence of the approach is hereby to choose as ba¬ 
sic unit cell of the lattice employed in the kMC 
simulations a suitable, large enough commensu¬ 
rate cell that embraces all primitive unit cells of the 
different lattices to be treated. This commensurate 
cell (and the lattice built up from it) then contains 
as sites simultaneously all sites of all lattices. 

In order to efficiently distinguish (during a sim¬ 
ulation) in which lattice state a given local region 
is currently in, we introduce an additional artifi¬ 
cial species (canonically named null) to repre¬ 
sent "disabled" sites. If say, a certain region is 
currently in one lattice structure, then the null 
species are placed on all sites corresponding to 
the other lattices in this region. A site "occupied" 
with a null species is hereby to be distinguished 
from an empty site. No elementary processes are 
available at sites occupied by a null species, i.e. 


elementary processes that require an empty site 
explicitly look for the sites to be empty and not 
"occupied" by null. In a corresponding region 
where all sites belonging to the other lattices are 
disabled through null species, multi-lattice Ip- 
kMC thus reduces to a regular lattice Ip-kMC sim¬ 
ulation and only the elementary processes corre¬ 
sponding to the Ip-kMC model of this lattice can 
possibly occur. 

In this setup, local transitions between the differ¬ 
ent lattices are finally enabled by introducing new 
elementary processes which involve the removal 
or placement of the null species at sites belong¬ 
ing to the lattices involved. Trivially, for the ox¬ 
ide reduction example this could be as simple as 
an elementary process placing the null species 
on one or more sites of the oxide lattice and re¬ 
moving the null species from nearby sites of the 
metal lattice. As shown below for the Pd oxide re¬ 
duction, local lattice transitions are unlikely that 
simple, but follow complex, multi-stepped reac¬ 
tion paths. The additional challenge of setting up 
and performing a multi-lattice Ip-kMC simulation 
as compared to traditional lattice Ip-kMC is there¬ 
fore to establish such atomistic pathways for the 
lattice transitions and implement them in form of 
corresponding elementary processes. 

From a purely algorithmic point of view as 
implemented in the kMC modeling framework 
kmos multi-lattice kMC is identical to single¬ 
lattice kMC which underlines that multi-lattice 
kMC is highly efficient. The only difference is that 
those elementary steps accounting for morpholog¬ 
ical changes typically involve significantly more 
sites (10-20 sites) than typical ordinary surface re¬ 
action steps (1-8 sites). The scaling with the num¬ 
ber of conditions (i.e. here sites) has been studied 
in a previous paper.- This shows that the computa¬ 
tional scaling with the number of conditions is sig¬ 
nificant, {ff{n\og{n) to 0’{n^)). However in terms 
of overhead there is no fundamental difference be¬ 
tween single-lattice kMC and multi-lattice kMC. 
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2.3 Literature Ip-kMC models of CO 
oxidation at the a/ 5-surface oxide 
and at Pd(lOO) 

Our multi-lattice Ip-kMC simulations describing 
the reduetion of the Pd(lOO) surfaee initially cov¬ 
ered with a \/5-surfaee oxide film build on existing 
Ip-kMC models for CO oxidation at the \/5-oxide 
and for CO oxidation at Pd(lOO). These single¬ 
lattice models and the elementary processes con¬ 
sidered in them have been established by extensive 
DFT calculations and have been deseribed in de¬ 
tail before.—!^ Here, we only briefly recapitulate 
their essentials for completeness and as a means to 
introduee the two lattice types. 



Figure 1: Illustration of the lattiee employed in 
the Ip-kMC model of CO oxidation at Pd(lOO). 
Shown is a top view of the Pd(lOO) surface with 
a primitive unit eell and its hollow (yellow eirele) 
and bridge (green rectangle) sites. CO molecules 
adsorb at bridge sites and oxygen molecules ad¬ 
sorb (dissoeiatively) at hollow sites. Here and 
henceforth Pd atoms in fce(lOO) positions are de¬ 
picted as large blue spheres, O atoms as small red 
spheres, and C as small grey spheres. 

Systematic DFT calculations identified the hol¬ 
low and bridge sites offered by the Pd(lOO) sur¬ 
faee, ef. Fig. [H as stable high-symmetry adsorp¬ 
tion sites for O and CO, respeetively.— The Ip- 
kMC CO oxidation model on Pd(lOO) eorrespond- 
ingly eonsiders all non-concerted adsorption, des¬ 
orption, diffusion, and Langmuir-Hinshelwood re- 
aetion processes involving these sites.— Oxygen 
molecules adsorb dissoeiatively on two empty 
next-nearest neighbor hollow sites,— CO ad¬ 
sorbs on an empty bridge site. Desorption pro- 


eesses are time-reversals of these non-aetivated 
adsorption proeesses, with rate constants ful¬ 
filling detailed balanee. Oxygen diffusion oe- 
curs between nearest-neighbor hollow sites, CO 
diffusion between nearest-neighbor bridge sites. 
Strong short-range repulsive interaetions between 
adsorbed reaetion intermediates as calculated by 
DFTi^iM are accounted for through site-bloeking 
rules that prevent processes leading to 0-0 pairs 
at nearest-neighbor hollow-hollow distanees, to 
CO-CO pairs closer or at next-nearest-neighbor 
bridge-bridge distance, and O-CO pairs at nearest- 
neighbor hollow-bridge distance.— CO oxidation 
proceeds via a Langmuir-Hinshelwood meeha- 
nism and at the temperatures of interest leads to an 
instantaneously desorbing CO 2 molecule. The rate 
constants for all these elementary processes have 
been computed with DFT-TST, specifically em¬ 
ploying the PBE exehange-eorrelation (xc) fune- 
tional.— Typical for CO oxidation at transition 
metal surfaees, O and CO diffusion barriers are 
significantly lower than the barriers of all other el¬ 
ementary proeesses involved. In the Ip-kMC sim¬ 
ulations this leads to frequent exeeutions of diffu¬ 
sion events at minute time increments. To increase 
the numerical efficiency and be able to simulate 
the extended time scales necessary to follow oxide 
reduction we artificially increased diffusion barri¬ 
ers for CO of 0.14eV and for O of 0.28 eV— by 
0.5 eV. This is permissible, if diffusion is then still 
fast enough to achieve an equilibration of the ad- 
layer ordering in between the other (rare) elemen¬ 
tary processes.^ We validated that this is indeed 
the ease, by performing the Ip-kMC simulations 
with diffusion barriers increased by 0.6 eV and 
0.7 eV without obtaining any significant changes 
in the oxide reduction rates reported below. 

As established by detailed first-principles statis¬ 
tical mechanics simulations by Rogal et al— only 
two of the high-symmetry adsorption sites offered 
by the intact \/5-oxide are relevant for the re¬ 
ducing environments considered here: a bridge 
and a hollow site as illustrated in Fig. [2l Equiv¬ 
alent to the Pd(lOO) ease, the Ip-kMC CO ox¬ 
idation model at the \/5-oxide correspondingly 
considers all non-concerted adsorption, desorp¬ 
tion, diffusion, and Eangmuir-Hinshelwood re¬ 
action processes involving these sites Eley- 
Rideal reaction processes as recently suggested by 
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Figure 2: Illustration of the lattice employed 
in the Ip-kMC model of CO oxidation at the 
\/5-oxide. Shown is a top view of the (\/5 x 
\/5)i?27° PdO(lOl) surface oxide reconstruction 
with a primitive unit cell of the oxide layer and 
its hollow (yellow circle) and bridge (green rect¬ 
angle) sites. O and CO adsorb at both site types. 
In the stoichiometric termination shown, O atoms 
occupy all hollow sites. Here and henceforth Pd 
atoms in surface oxide positions are depicted as 
large green spheres. 


Hirvi et ah— are in principle also considered, but 
do not play a role at the low CO pressures ad¬ 
dressed here. Adsorption is non-activated, apart 
from dissociative adsorption of O 2 into neighbor¬ 
ing bridge sites, which is hindered by a sizable ad¬ 
sorption barrier of 1.9 eV. As for the Pd(lOO) Ip- 
kJVlC model, all Ip rate constants of this \/5-oxide 
model have been computed with DFT-TST, em¬ 
ploying the PBE xc functional.— DFT-computed 
nearest-neighbor lateral interactions are taken into 
account in these rate constants as detailed by Ro- 
gal et al. 

2.4 Computational DFT and kMC 
settings 

The multi-lattice Ip-kMC simulations of the ox¬ 
ide reduction process consider additional ele¬ 
mentary processes to model the atomistic path¬ 
way from oxide to metal lattice. This pathway 
and the involved elementary processes are de¬ 
scribed in Sections [T2ll3. 41 below. Their Ip rate 
constants have been calculated with DFT em¬ 
ploying the same TST-based expressions— and 


the same PBE xc functional— as used in the 
Pd(lOO) and \/5-oxide Ip-kMC models. Employ¬ 
ing the same computational setup as detailed in 
ref. I 28 L the DFT calculations have been performed 
with the CASTEP package— via the Atom¬ 
ic Simulation Environment.— Previous 
work— established that the employed standard li¬ 
brary ultrasoft pseudopotentials accurately repro¬ 
duce the full-potential \/5-oxide adsorption ener¬ 
getics of Rogal et al. Geometry optimizations 
and transition-state searches were done with super¬ 
cell geometries containing asymmetric slabs with 
one (frozen) layer of Pd(lOO) and one layer of 
PdO(lOl), separated by at least 10 A vacuum. Ge¬ 
ometry optimizations were performed using the 
BFGS method with a force threshold of 0.05 eV/A. 
The energetics of the thus obtained intermediate 
geometries were refined in single-point calcula¬ 
tions with further Pd(lOO) layers added to the 
slab. Transition state searches have been per¬ 
formed with the climbing-image nudged elastic 
band method— or simple drag procedures with the 
nature of the transition state confirmed through 
frequency calculations. From test calculations 
with additional Pd(lOO) layers added we estimate 
the uncertainty in the thus determined barrier val¬ 
ues as 0.1 eV. The gas-phase calculations for O 2 
and CO were performed in a cubic supercell of 
15 A box length and F-point sampling. At a plane- 
wave cutoff of 400 eV and a [4 x 4 x 1] k-point 
grid for the commensurate surface oxide unit cell 
and a [2 X 2 X 1] k-point grid for the (2 x 2) sur¬ 
face oxide unit cell in the slab calculations conver¬ 
gence tests demonstrate a numerical convergence 
of the obtained DFT barriers and binding energies 
to within 0.1 eV. 

All Ip-kMC simulations have been performed 
using the kmo s framework.- The periodic bound¬ 
ary simulation cell contains (40 x 20) \/5-unit 
cells and thus 1600 sites in case of the single¬ 
lattice Ip-kMC simulations on the surface oxide. 
Systematic tests with varying simulation cell sizes 
showed the obtained oxide reduction rates to be 
fully converged at this setting. In the multi-lattice 
case one commensurate unit-cell contains already 
a total of 19 sites from both the metal and the 
oxide lattice. We employed periodic boundary 
simulation cells comprising (10 x 20 ), (20 x 20 ) 
and (40 x 20) such commensurate unit-cells to 
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specifically assess the role of elementary processes 
across oxide and metal domains (see below). In all 
cases, averages over 10 kMC trajectories starting 
with different random number seeds were taken to 
obtain the correct transient quantities. 

3 Results 

3.1 Single-lattice Ip-kMC simulations 
of oxide reduction 

As mentioned initially the central objective of 
our simulations is to analyze the oxide reduction 
rates obtained in recent high-resolution X-ray pho¬ 
toelectron spectroscopy experiments by Fernan¬ 
des et ah— Within a small range of temperatures, 
they reported large variations for the total time re¬ 
quired to fully reduce an initially prepared \/5- 
oxide layer by exposure to a pure CO atmosphere 
of 5 X 10^^^bar. Specifically, at 303K the ox¬ 
ide was reduced after approximately 90 minutes, 
while the same process took only approximately 
12 min at 343 K and less than 4 min at 393 K. Pre¬ 
vious theoretical work addressing the surface ox¬ 
ide stability employed the Ip-kMC approach with 
a single lattice, namely the \/5-oxide Ip-kMC 
model. They monitored the coverage of O in the 
hollow sites as an indicator for the oxide stabil¬ 
ity.^^ Using this criterion Ip-kMC simulations fail 
qualitatively to reproduce the experimental find¬ 
ings of Fernandes et al.. Even at the highest tem¬ 
perature of 393 K no reduction of the oxide is ob¬ 
served at all after 120 mins, i.e. the oxygen cover¬ 
age remains essentially intact over the entire time. 
This can be explained by the lack of CO on the 
surface, which implies that an efficient removal of 
oxygen via CO oxidation is not possible. 

Since an error in the reported experimental con¬ 
ditions large enough to explain this discrepancy 
seems unlikely we focus on the uncertainty in the 
DFT parameters deployed in the model. A pos¬ 
sible reason for this discrepancy could be that 
the relevant CO adsorption energy of —0.92eV 
as calculated by DFT-PBE is too small, and/or 
the CO reaction barrier towards CO 2 formation of 
0.9 eV^^ is too high. We explore these possibilities 
by running the Ip-kMC simulations using modi¬ 
fied rate constants. Specifically, we vary the CO 



Eigure 3: Single-lattice kMC simulations for 
\/5-oxide reduction at the three temperatures em¬ 
ployed in the experiments and using the O cov¬ 
erage at the hollow sites (0Oh) indicator for 
the progressing reduction of the surface oxide (see 
text). The Ip rate constants are modified to achieve 
at least some agreement with the experimental 
findings: (a) CO adsorption energy lowered by 
0.2 eV. The reduction times for all three tempera¬ 
tures are in the right order of magnitude, but the or¬ 
dering with temperature is wrong, (b) CO adsorp¬ 
tion energy lowered by 0.5 eV; CO oxidation bar¬ 
rier increased by 0.15 eV. The oxide reduction time 
correctly decreases with temperature, but the re¬ 
duction time at 303 K disagrees qualitatively with 
experiment. 
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oxidation barrier and the CO adsorption energy on 
the bridge and hollow sites within ±leV. How¬ 
ever, no sueh kMC simulation reproduees the ex¬ 
perimentally observed behavior: Either the time 
seale for reduetion is too long or too short for all 
three temperatures, or if the simulations exhibit a 
eorreet time seale for at least one of the temper¬ 
atures, then the trend with temperature is wrong, 
i.e. the reduetion time does not deerease for eaeh 
higher temperature as found in the measurements. 
We illustrate this in Fig. [3] for two oases where 
at least partial agreement with experiment is ob¬ 
tained: When the CO adsorption energy is low¬ 
ered by 0.2 eV and the CO reaetion barrier is left 
unohanged, reasonable reduetion time soales are 
obtained for all three temperatures, but the order¬ 
ing is wrong: At the intermediate temperature of 
343 K reduetion is slower than at the lower temper¬ 
ature of 303 K. When instead the CO adsorption 
energy is lowered by 0.5 eV and the CO oxidation 
barrier is simultaneously inereased by 0.15 eV, this 
ordering is eorreet, but the time seale for oxide re¬ 
duetion at 303 K is mueh too long. 

It is diffieult to generally rule out that also vary¬ 
ing the rate eonstants of all other elementary pro- 
eesses would not eoineidentally yield a parame¬ 
ter set that reproduees the experimental reduetion 
times. Nevertheless, the failure to obtain sueh a 
parameter set when modifying the two proeesses 
that most obviously relate to the reduetion time, let 
us believe that a shorteoming of the employed Ip 
rate eonstants is unlikely the reason for the quali¬ 
tative diserepaney to the measurements. We there¬ 
fore tend to eonelude that the reduetion of the y/S- 
oxide layer does not happen at intaet oxide patehes 
(for whieh the O eoverage at the hollow sites is a 
useful stability indieator). 

3.2 Atomistic pathway for the ^/ 5 - 
oxide reduction 

In the eontext of oxide formation and reduetion 
effieient eross-reaetions at the boundary between 
oxide and metal domains have long been dis- 
eussed.— Fernandes et al. reported long induetion 
times of the order of minutes at all measured tem¬ 
peratures.— This renders a erueial role of eross- 
reaetions at extended surfaee imperfeetions like 
steps or disloeations rather unlikely. Instead we 



Figure 4: Top view of the supereell model em¬ 
ployed in the oxide reduetion DFT ealeulations, 
spanning four \/5-oxide unit eells in a (2 x 2) ar¬ 
rangement. Eaeh intaet \/5-oxide unit eell eon- 
tains four O atoms, labeled from A to D, and four 
Pd atoms, two of whieh (marked through hateh- 
ing) are loeated in hollow sites of the underlying 
Pd(lOO) substrate. The arrows in the upper right 
three unit-eells indieate the direetions in whieh the 
two Pd atoms ean be shifted to generate a minimal 
Pd(lOO) pateh: This defines the three eases I (up¬ 
per left unit-eell), II (upper right unit-eell) and III 
(lower right unit-eell) referred to in the text. 
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suspect that a slow initial reduction of ideal ox¬ 
ide terraces locally generates minimal Pd(lOO) 
patches that then enable an efficient continued re¬ 
duction through cross-reactions. Admittedly, also 
for this suggestion the observed induction time is 
somewhat long. In order to scrutinize this picture 
with multi-lattice Ip-kMC simulations we there¬ 
fore establish an atomistic pathway from the in¬ 
tact \/5-oxide to reduced Pd(lOO) and focus on 
the structural evolution of the surface oxide film 
when progressively removing O atoms. These 
studies are consistently performed with large su¬ 
percell models spanning four \/5-oxide surface 
unit cells in a (2 x 2) arrangement as shown in 
Fig. m This model is large enough to study the 
effects of O removal on one y/S cell that is still 
surrounded by intact \/5-oxide units in every di¬ 
rection. For these studies, an important structural 
aspect of the commensurate \/5-oxide/Pd(100) in¬ 
terface is that in every \/5-oxide unit cell two of 
the four Pd atoms in the oxide layer are located 
in fourfold hollow sites of the underlying Pd(lOO) 
substrate, with the other two being rather close to 
bridge sites, cf. Fig. |4l After complete reduction, 
the Pd atoms of the surface oxide will form do¬ 
mains of a new Pd(lOO) layer,— in which they 
are also located in the fourfold hollow sites. This 
means that during the local oxide decomposition 
two out of four Pd atoms in the oxide layer do not 
necessarily have to significantly change their lat¬ 
eral position. 

In the stoichiometric \/5-oxide termination there 
are four O atoms per \/5 unit cell. Two are located 
above the Pd atoms of the oxide layer (occupying 
the hollow sites described above) and two below, 
as shown in Fig. HI In agreement to previous DFT 
studies performed in smaller supercells,- we find 
no significant structural changes upon removal of 
any single of these four O atoms in one y/S unit 
cell of our model. In particular, none of the oxide 
Pd atoms change their lateral registry. The same is 
observed, when removing pairs of O atoms in all 
possible combinations. Apparently neither single 
O vacancies, nor O divacancies lead to an unacti¬ 
vated destabilization of the \/5-oxide structure. 

In this situation, we focus on a straightforward 
way to generate a minimal Pd(lOO) patch within 
the surface oxide layer. We explore if this (to¬ 
gether with an increasing number of O vacancies) 


does not constitute a pathway with at least a very 
low activation barrier. This pathway centers on 
the above mentioned realization that two of the 
four Pd atoms in the \/5-oxide unit cell are al¬ 
ready located in fourfold hollow sites with respect 
to the underlying Pd(lOO) substrate. A direct way 
to convert one entire y/5 cell into a Pd( 100)-type 
patch would therefore be to simply shift the other 
two Pd atoms from their close to bridge positions 
to adjacent hollow sites. As illustrated in Fig. |4] 
there are three ways to realize this: Either both 
Pd atoms move to hollow sites in the same direc¬ 
tion (cases I and II), or the two Pd atoms move 
to hollow sites in opposite directions (case III). If 
we enforce such lateral shifts on the Pd atoms of 
one unit cell for the intact surface oxide, they re¬ 
lax back to their original positions upon geometry 
optimization. The same holds, if this is done with 
one O vacancy in the \/5-oxide unit cell, or when 
moving only one of the two Pd atoms. 

In contrast, metastable structures result when the 
same movements are applied in the presence of an 
O divacancy. Under reaction conditions such an O 
divacancy could arise from the random reduction 
of neighboring O atoms in the oxide but possibly 
also due to a consequential reaction step next to 
a metal to surface oxide boundary. Following a 
systematic development from the fully intact sur¬ 
face oxide similar to a cluster expansion divacan¬ 
cies are an obvious next candidate motif after sin¬ 
gle vacancies. Specifically this is when removing 
the two upper O atoms, labeled A and B in Fig. HI 
that are directly coordinated to one of the moving 
Pd atoms. For the three cases defined in Fig. HI 
we calculate the relative energies of the resulting 
metastable structures compared to the original di¬ 
vacancy structure with the Pd atoms in their regu¬ 
lar positions. They are —0.40eV (I), —0.25 eV (II) 
and -hO. 17 eV (III). Cases I and II lead therefore to 
more stable intermediate structures. In particular 
for the most stable case I we furthermore find only 
a negligibly small activation barrier below 0.1 eV. 

This suggests the following initial steps in the 
oxide reduction pathway: Upon local formation 
of an O divacancy as shown in Fig. [5^, the two 
bridge-site Pd atoms of the corresponding \/5- 
oxide cell shift laterally and essentially unacti¬ 
vated to yield the intermediate structure shown in 
Fig. [ 5 J 5 . This metastable structure is more stable 
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Figure 5: Top view of intermediate struetures along the suggested oxide reduetion pathway, (a) Divaeaney 
ereated by removal of two neighboring upper O atoms in one \/5-oxide unit-eell (labeled A and B in Fig. 
IH). (b) Intermediate strueture formed from (a) through an essentially unaetivated shift of the two bridge- 
site Pd atoms as in ease I of Fig. HI This strueture is by 0.4 eV more stable than strueture (a), (e) Pd(lOO) 
pateh formed from (b) after the aetivated diffusion of the sub-surfaee O labeled D in Fig. |4]to the fourfold 
on-surfaee site in the eenter of the thus formed Pd(lOO) pateh (see text). 


than the divaeaney by 0.4 eV. It already features a 
minimal Pd(lOO) pateh with four Pd atoms in hol¬ 
low sites, but still eontains the two O atoms of the 
\/5-oxide eell that were originally below the Pd 
layer. One of these two, labeled C in Fig. HI auto- 
matieally moves to a position above the Pd atoms 
of the Pd(lOO) pateh along with the lateral shift 
of the bridge-site Pd atoms. The other O atom 
D remains in a sub-surface position as depicted in 
Fig.[5j). 

We calculate a barrier of 0.66 eV for the diffu¬ 
sion process that lets this O atom pop up to the on- 
surface site E in the center of the Pd(lOO) patch, 
cf. Fig. [St. There the O atom is by 0.11 eV more 
stable than in its original sub-surface site. From 
there and as further detailed below, the on-surface 
O atom can be reacted off through Langmuir- 
Hinshelwood-type CO oxidation. Alternatively, it 
can diffuse away. Both cases leave a bare mini¬ 
mal Pd(lOO) patch as the end product of the here 
suggested oxide reduction pathway. Continued re¬ 
duction of the oxide layer will then lead to an in¬ 
creasing number of such Pd(lOO) patches. As the 
Pd density in the intact surface oxide is 20 % lower 
than in a Pd(lOO) layer, Pd vacancies will emerge 
along with the formation of these patches. The 
diffusion of such Pd vacancies will be very fast. 
The calculated activation barrier of only ~ 0.2eV 
is similar to the diffusion of Pd on Pd(lOO).— 
Quickly, this fast diffusion will lead to the coales¬ 
cence of the individual Pd(lOO) patches and with 


that to the formation of extended Pd(lOO) islands 
as seen in experiment.— 

3.3 Elementary processes at the 
Pd(lOO) patch under reduction 
conditions 



Figure 6: Top view of the Pd( 100) patch as formed 
during the oxide reduction process. Shown are all 
high-symmetry binding sites explored to assess the 
binding properties of the patch itself (sites A-E) 
and at its boundary (sites E-K). Pd atoms of the 
formed patch are darkened. 

The Pd(lOO) patch is the central new feature 
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within the obtained oxide reduetion pathway. Dur¬ 
ing initial oxide reduction it may emerge upon suf¬ 
ficient local oxygen depletion of the \/5-oxide. Its 
existence may then enable new elementary pro¬ 
cesses specific to either the patch itself or the 
boundary between patch and surrounding oxide. 
As the next step we therefore proceed to a sys¬ 
tematic identification and description of corre¬ 
sponding processes. As starting point for this we 
first consider the binding at the different high- 
symmetry sites illustrated in Fig. 0 This com¬ 
prises sites on top of the patch and at its bound¬ 
ary. Specifically, sites labeled A-D and site E cor¬ 
respond to Pd(100)-type bridge and hollow sites 
on the patch, respectively. Sites F-J are sites di¬ 
rectly at the patch/oxide boundary, whereas site K 
corresponds to a \/5-oxide bridge site in the im¬ 
mediate vicinity of the patch. Table [U summarizes 
the computed binding energies at these sites for the 
two "oxidation" stages of the patch along the oxide 
reduction pathway: With the sub-surface oxygen 
atom still beneath it, cf. Fig. [5j), and without it. 
The latter situation results as a consequence of the 
activated diffusion of the O atom to the on-surface 
hollow site of the patch and subsequent reaction or 
diffusion, vide infra, and is the situation depicted 
in Fig. [6l 

With the exception of site B directly adjacent 
to the sub-surface O atom, the effect of the latter 
atom on the binding properties of the sites on top 
of the patch is rather small. Quite similar bind¬ 
ing energies are obtained for these sites for the 
two cases. Cum grano salis these binding ener¬ 
gies on the patch sites are also quite similar to 
those computed for the extended Pd(lOO) surface. 
Roughly we obtain a preferred O binding at the 
hollow site E of about — 1.25eV. This should be 
compared to — 1.25 eV at the corresponding site on 
Pd(lOO). In turn, the calculated preferred CO bind¬ 
ing at the bridge sites A-D of about ~ — 1.4-2.1 eV 
should be compared to —1.93 eV at bridge sites of 
Pd(lOO). Without the presence of the sub-surface 
O atom this also extends to the bridge sites F and 
H at the interface between the patch and the \/5- 
oxide. Again the CO binding properties are very 
similar to regular bridge sites at extended Pd(lOO). 
Furthermore considering the close geometric and 
energetic vicinity of these different sites we there¬ 
fore also expect similarly low diffusion barriers as 


Table 1: Calculated DFT binding energies for 
CO and O at high-symmetry sites on top of the 
Pd(lOO) patch and at its boundary. Shown are re¬ 
sults for the situation with the sub-surface O atom 
still beneath the patch (labeled "w/ sub-surf O") 
and without the sub-surface O atom (labeled "w/o 
sub-surf O"). Sites A-K are defined in Fig. 0 
Fields marked with a letter indicate that the adsor¬ 
bate moved to the respective site during geometry 
optimization. Fields marked " indicate sites for 
which no binding was obtained. Binding in sites 
marked with an asterisk ended up in a threefold 
coordination after geometry optimization. All en¬ 
ergies are in eV. 



w/ sub-surf 0 

w/o sub-surf 0 


CO 

0 

CO 

0 

A 

-1.72 

E 

-1.85 

E 

B 

— 

— 

-1.73 

I 

C 

-1.40 

D 

-1.62 

-0.98* 

D 

-2.06 

-0.96 

-2.06 

-1.00 

E 

D 

-1.22 

-1.90* 

-1.25 

E 

-1.49 

-0.96 

-1.72 

-1.25 

G 

-1.27 

-0.40 

-1.31 

-0.47 

H 

-0.83 

J 

-1.56 

J 

I 

— 

— 

B 

-0.96 

J 

H 

-0.39 

H 

-0.75 

K 

-1.03 

-0.30 

-1.10 

-0.37 
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on Pd(lOO). Correspondingly, there will be a pre¬ 
dominant presenee of the adsorbates in the most 
favorable sites, i.e. CO in the bridge site D and O 
in the hollow site E at the temperatures employed 
in the oxide reduetion experiments. 

Multiple adsorption on one pateh in form of CO 
adsorption at several bridge sites, or O and CO 
eoadsorption on hollow and bridge site, respee- 
tively, is not possible aeeording to our DFT eal- 
eulations. Also in this respeet, adsorption at the 
Pd(lOO) pateh seems very similar to adsorption 
at extended Pd(lOO). Likewise the disturbanee of 
the \/5-oxide ereated by the Pd(lOO) pateh seems 
quite short-ranged. Already at the oxide bridge 
site K direetly adjaeent to the pateh, O and CO 
binding is essentially identieal to the binding at the 
bridge site on the stoiehiometrie y/S termination. 
The CO binding is ealeulated as — 1.03eV and 

— l.lOeV at site K (with and without sub-surfaee 
O underneath the adjaeent Pd(lOO) pateh) and as 

— 1.08 eV on the intaet \/5-oxide. The O binding 
energy is 0.30 eV, —0.37 eV and —0.19eV for the 
three eases, respeetively. Finally, the binding at the 
interfaeial sites F, G, H and I in the presenee of the 
sub-surfaee oxygen atom and at the interfaeial site 
H without the sub-surfaee O atom is intermediate 
to the one obtained on the pateh and on the intaet 
\/5-oxide. 

In light of these results, we ean now systemat- 
ieally analyze the new elementary proeesses that 
may arise due to the pateh under the redueing eon- 
ditions of the experiment and assess their qualita¬ 
tive relevanee for the oxide reduetion proeess. CO 
adsorption and desorption proeesses ean take plaee 
at the pateh bridge sites, with a maximum eover- 
age of one CO moleeule on a pateh. This is re¬ 
gardless of whether the sub-surfaee O atom is still 
present underneath the pateh or not. CO diffusion 
on the pateh itself will be fast, and will lead to a 
preferential population of bridge site D at the tem¬ 
peratures employed in the reduetion experiments. 
After eoaleseenee of nuelei to larger Pd(lOO) do¬ 
mains, equivalent proeesses (CO adsorption, des¬ 
orption, diffusion) will oeeur, with properties as 
on extended Pd(lOO). In prineiple, CO diffusion 
eould also oeeur from pateh (or Pd(lOO) domain) 
to oxide or viee versa. An important point to real¬ 
ize here, however, is that CO binding at the pateh 
(and at Pd(lOO)) is about 1 eV more stable than at 


the hollow or bridge sites of the \/5-oxide. This al¬ 
ready imposes a quite large mere thermodynamie 
barrier for CO diffusion onto the oxide. Simulta¬ 
neously CO diffusion proeesses from the oxide to 
the pateh will be rare, owing to the near zero CO 
eoverage on the surfaee oxide under the eonditions 
of the reduetion experiments, ef. Seetion lXTl 

If all sites A-E of the pateh are unoeeupied, the 
sub-surfaee O atom ean pop up into the on-surfaee 
fourfold hollow site E of the pateh, with a eal¬ 
eulated barrier of 0.66 eV. The reverse proeess is 
possible with a eorresponding barrier of 0.77 eV. 
From the pateh site E the oxygen atom ean dif¬ 
fuse to an adjaeent vaeant hollow site of the sur¬ 
rounding oxide, or viee versa. Sinee O binding is 
about 0.4 eV more stable at the oxide hollow sites, 
this diffusion will preferentially take oxygen away 
from the pateh (or eoaleseed larger Pd(lOO) do¬ 
mains). For O diffusion between the oxide hollow 
sites direetly adjaeent to the pateh (those above site 
K in Fig. we find a slightly lowered diffusion 
barrier of 1.1 eV as eompared to the 1.4 eV for the 
diffusion over the intaet surfaee oxide. This lower¬ 
ing ean be rationalized with the higher geometrie 
flexibility of the underlying Pd atoms in the par¬ 
tially redueed surfaee oxide environment. Oxy¬ 
gen diffusion from site E to adjaeent bridge sites 
of the \/5-oxide faee a similarly large thermody¬ 
namie barrier as the CO diffusion proeesses due to 
the mueh weaker O binding at these oxide sites. 
As these bridge sites are furthermore essentially 
all unoeeupied under the eonditions of the redue¬ 
tion experiments, the reverse diffusion proeesses 
bringing O atoms from adjaeent oxide bridge sites 
to the pateh will also oeeur at a negligibly low 
rate. Assoeiative oxygen desorption requires two 
nearby O atoms. It ean thus only take plaee either 
after eoaleseenee of nuelei (whieh ereates next- 
nerest neighbor Pd(100)-type hollow sites) or in 
ease of an individual pateh by involving one O 
atom from the pateh site E and one O atom from 
nearest neighbor O hollow sites of the oxide. Due 
to the similar binding energies of O on both metal 
and oxide, sueh proeesses are as rare as they were 
in the single-lattiee Ip-kMC simulations foeusing 
only on the surfaee oxide. Correspondingly, sueh 
proeesses are unlikely to have a relevant effeet on 
the reduetion ki netics. 

Finally, new CO oxidation reactions can be en- 
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Figure 7: Top view of the Pd(lOO) pateh illus¬ 
trating possible CO oxidation reaetions aeross the 
pateh/oxide boundary. The CO moleeule is ad¬ 
sorbed at the most stable bridge site D of the 
Pd(100)-patch. The calculated barriers for the four 
reaction processes are: NW 1.63 eV, SW 1.46 eV, 
SE 1.54eV, and E 0.95 eV. 


abled at the patch or its boundary. One possibil¬ 
ity is a reaction of a CO molecule adsorbed on 
the patch with the sub-surface O atom. Eor a CO 
molecule occupying the preferred bridge site D we 
calculate a barrier of 1.27 eV for this process. This 
is rather high for this process to play a significant 
role at the temperatures employed in the experi¬ 
ments. Alternatively, the CO molecule could react 
with O atoms of the surrounding \/5-oxide. Eo- 
cusing on the preferred binding site D Eig. |7] illus¬ 
trates the four corresponding reactions to the clos¬ 
est such O atoms and compiles the calculated bar¬ 
riers for these processes. Only one of these pro¬ 
cesses exhibits a barrier below 1 eV and should 
therewith exhibit a similar reactivity as the reac¬ 
tion processes occurring on either Pd(lOO) or on 
the \/5-oxide. As last possibility, an O atom oc¬ 
cupying the on-surface site E of the patch could 
react with a nearby adsorbed CO molecule. Due 
to the near zero CO coverage on the oxide during 
the reduction experiments, such a process is un¬ 
likely to occur with a CO molecule adsorbed on 
the y/S though. At the same time, the strong re¬ 
pulsive interactions prevent CO adsorption at the 
bridge sites A-D of the patch, if an O atom occu¬ 


pies the hollow site E. This leaves as only possi¬ 
bility for such a cross-reaction, the reaction with a 
CO molecule adsorbed at the bridge site E, cf. Eig. 
|6l In general, such a reaction process can not have 
a strong influence on the oxide reduction process 
though. As O diffusion from the \/5-oxide onto 
the patch is rare, this process will predominantly 
occur after the sub-surface O atom has popped up 
to the surface. The process then takes place once 
per formed patch and is thereafter unlikely to re¬ 
occur. 

Summarizing these insights, the qualitatively 
new feature emerging from the presence of a 
formed Pd(lOO) patch during the oxide reduction 
process is the significantly enhanced stabilization 
of CO at the patch as compared to the \/5-oxide. 
Whereas the CO coverage on the latter is close to 
zero under the reducing conditions of the experi¬ 
ment, a much higher CO coverage can thus be ex¬ 
pected at formed patches. An oxidation reaction 
of such CO molecules across the patch/oxide in¬ 
terface with O atoms of the surface oxide then pro¬ 
vides a much more efficient means to deplete the 
surface oxygen compared to the two mechanisms 
possible at the \/5-oxide alone. There the asso¬ 
ciative O desorption is a rare process due to the 
strong O binding energy on the oxide. Similarly 
the CO oxidation is a rare process due to the low 
CO coverage at the surface. In contrast, the new 
reduction process enabled at the oxide/patch inter¬ 
face only depends on the CO adsorption, which 
is comparatively frequent in the reducing environ¬ 
ment of the experiment, and on the replenishment 
of the created O vacancies at the oxide sites at the 
interface. The latter occurs via O diffusion on the 
surface oxide, which is not a very fast process. Yet, 
it is still faster than the oxide-only reaction mech¬ 
anism, considering in particular that the number 
of formed Pd(lOO) patches will scale with the in¬ 
creasing number of concomitantly created O va¬ 
cancies in the oxide layer. 

3.4 Multi-lattice Ip-kMC model 

We now proceed by casting the established pic¬ 
ture of the reduction process via the formation of 
Pd(lOO) patches into an explicit multi-lattice Ip- 
kMC model. Eor regions of Pd(lOO) and intact 
\/5-oxide this multi-lattice Ip-kMC model will 
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employ the literature single-lattice Ip-kMC mod¬ 
els as before. The new elementary processes aris¬ 
ing in the context of the Pd(lOO) patches and their 
DFT calculated barriers are listed as follows: 

A local reduction begins whenever an oxygen di¬ 
vacancy involving two neighboring hollow sites on 
the \/5-oxide, cf. Fig. [5l is formed during the 
Ip-kMC simulation. The patch formation is then 
modeled as a non-activated process, i.e. it auto¬ 
matically follows the formation of the divacancy. 
In terms of the multi-lattice representation the 
patch formation corresponds to deactivating the 
two bridge and two hollow sites associated with 
the corresponding \/5-oxide unit cell by changing 
their occupation status to null. In their place, 
Pd(lOO) lattice sites are activated. Specifically, 
this concerns Pd(lOO) bridge sites at sites labeled 
A, B, C, D, F, H in Fig. the occupation of which 
is changed from null to empty. Furthermore, 
a Pd(lOO) hollow site is created empty at the 
site labeled E in Fig. 0 The quenched-in oxygen 
atom is modeled to remain at its sub-surface y/S- 
oxide site. However, now that the Pd(lOO) hollow 
site labeled E in Eig.[^is created, this sub-surface 
O atom can hop there with a barrier of 0.66 eV. 
However, this only, if the Pd(lOO) hollow site 
and its four nearest-neighbor Pd(lOO) bridge sites 
are empty in accordance with the site-blocking 
rules used in the Pd(lOO) single-lattice Ip-kMC 
model.— Erom there, the oxygen atom can reverse 
the diffusion with a barrier of 0.77 eV. Alterna¬ 
tively, if there are already other immediately adja¬ 
cent Pd(lOO) patches, it can also regularly diffuse 
to corresponding nearest-neighbor Pd(lOO) hollow 
sites. O diffusion from the patch to surrounding 
\/5-oxide (and the reverse back-diffusion process) 
exhibits a barrier of 0.5 eV (0.9 eV). The hollow- 
hollow diffusion of oxygen on the \/5-oxide im¬ 
mediately adjacent to a Pd(100)-patch has a bar¬ 
rier of 1.1 eV, and is therewith accelerated as com¬ 
pared to diffusion on the intact y/5 with its barrier 
of 1.4 eV. 

The second oxygen atom that is originally in the 
sub-surface \/5-oxide site and pops up to the sur¬ 
face during the formation of the Pd(lOO) patch, cf. 
Eig. ISh, has a binding energy of — 1.42eV, irre¬ 
spective of whether the other adjacent quenched- 
in sub-surface O atom is still present or not. Even 
though thus exhibiting a somewhat stronger bind¬ 


ing by 0.17 eV, we still approximate this as a reg¬ 
ular O atom in a Pd(lOO) hollow site. In the 
multi-lattice scheme we correspondingly activate 
another Pd(lOO) hollow lattice site at the location 
that is obtained by going from site E one lattice 
constant to the left in Eig. and directly occupy 
it with the popped-up O atom. Simultaneously, 
the sub-surface \/5-oxide site in which the oxygen 
atom resided before the patch formation is deacti¬ 
vated by setting its occupation to null. 

If the new patch lies next to already existing 
patches, CO diffusion from neighboring Pd(lOO) 
sites onto bridge sites of the new Pd(lOO) patch 
is possible, if the nearest-neighbor hollow as well 
as all up to next-nearest-neighbor bridge sites are 
empty - again in accordance with the site-blocking 
rules for CO diffusion in the Pd(lOO) single-lattice 
Ip-kMC model.— If all Pd(lOO) patch bridge and 
hollow sites are empty, non-activated CO adsorp¬ 
tion is possible onto a bridge site. CO desorption 
has a desorption barrier of 1.93 eV and obeys de¬ 
tailed balance with the adsorption process as in 
the Pd(lOO) single-lattice Ip-kMC model. As ad¬ 
ditional reaction processes we only consider the 
most likely process labeled E in Eig. |7]with a bar¬ 
rier of 0.95 eV, involving a CO molecule adsorbed 
at site D of the patch and an O atom adsorbed at a 
directly adjacent upper hollow \/5-oxide site. 

With these prescriptions and new processes the 
multi-lattice Ip-kMC simulation will allow for the 
generation of Pd(lOO) patches, which if generated 
next to each other will increasingly lead to larger 
Pd(100)-type domains in the course of a reduc¬ 
tion condition simulation. As a Pd(lOO) unit-cell 
of the commensurate \/5-oxide/Pd(100) interface 
exhibits only 4/5 of the Pd atom density, these 
patches will not correspond to dense Pd(lOO) ter¬ 
races though. Instead, they will contain a regu¬ 
lar array of Pd vacancies. In reality, fast Pd self¬ 
diffusion will lead to the coalescence into closely 
packed Pd(lOO) islands, and connected with it to 
an entire class of additional elementary processes. 
However, in the here proposed atomistic pathway, 
cf. Pig.m the Pd(100)-patch formation step shifts 
the two mobile Pd atoms in such a way that the 
emerging Pd(lOO) vacancy is located furthest from 
the relevant \/5-oxide/Pd(100) domain boundary, 
where the cross-reaction E of Eig. |7]can take place. 
As such we anticipate that the coalescence of Pd 
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vacancies through Pd self-diffusion oeeurs pre¬ 
dominantly displaced from the reduetion front and 
then has no effeet on the aetual reduetion kineties. 
Within this view, we neglect Pd self-diffusion in 
the multi-lattiee Ip-kMC model, sueh that the end 
produet of a eomplete reduction of the \/5-oxide 
film eorresponds to a Pd(lOO) surface with a regu¬ 
lar array of Pd vaeancies. 

Table 2: DFT parameters entering rate constants 
of elementary proeesses near Pd(100)/PdO(\/5 x 
\/5)R21° domain boundary. All energies are in eV. 


process 

barrier 

oxide reeonstruetion 

0.0 

0 spillover 

0.66 

0 reverse spillover 

0.77 

0 hollow diffusion near nucleus 

1.1 

eross-reaetion (E) 

0.95 


3.5 Multi-lattice Ip-kMC simulations 
of oxide reduction 

Within the multi-lattiee Ip-kMC model we now 
address the oxide reduction experiments of Fer¬ 
nandes et al. The qualitatively new feature in¬ 
troduced by the explieit aeeount of reduced ox¬ 
ide patehes in the multi-lattice simulations is the 
possibility of reduetion through eross-reactions at 
\/5-oxide/Pd(100) domain boundaries. Within the 
picture of some initial slow reduction during the 
long induction time seen in the experiments, we 
employ kMC simulation cells, in which initially 
one row of \/5-oxide along the [010]-direction is 
already removed, cf. Fig.[^ To assess the sensitiv¬ 
ity on this point, we conduct the multi-lattice Ip- 
kMC simulations using (10 x 20), (20 x 20) and 
(40 X 20) eells, thereby systematieally varying the 
width of the initially intaet \/5-oxide domain. As 
demonstrated in Fig. |9]we obtain only small vari¬ 
ations in the oxide reduetion times. These do not 
affeet at all the eonelusions put forward below and 
confirm the understanding that the oxide reduetion 
rates seen in experiment are governed by a fast re¬ 
duction following some slow aetivation during the 
induetion period. 

Figure |9] shows the surface oxide coverage as 
a function of time using the experimental eon- 



Figure 8: Illustration of the initial eonfiguration 
used in the multi-lattice Ip-kMC simulations. For 
clarity, only a small (4x8) commensurate \/5- 
oxide eell area of the larger simulation eell is 
shown. One row of \/5-oxide eells has been re- 
dueed to Pd(lOO) nuclei along the (100) direction 
of the multi-lattice kMC model. Shown in paren¬ 
thesis are the corresponding erystallographie axes 
of the PdO erystal. 
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Figure 9: Multi-lattice Ip-kMC simulations of the 
\/5-oxide reduction at the three experimentally 
employed temperatures. Shown is the fraction of 
the surface covered by the surface oxide as a func¬ 
tion of time using a (20 x 20 ) (solid, bold lines), 
a (10 X 20) (dashed, thin lines), and a (40 x 20) 
(dotted, thin lines) simulation cell. The slight dif¬ 
ference in the initial \/5-oxide fraction for the dif¬ 
ferent cell sizes results from the fact that in every 
cell initially one row of oxide is removed. 

ditions of 5 x 10 bar CO pressure, and tem¬ 
peratures of 303 K, 343 K, and 393 K. In con¬ 
trast to the single-lattice Ip-kMC simulations dis¬ 
cussed in Section 13.1[ very good agreement with 
the experimentally observed times and trends is 
obtained - without the need to modify any of the 
Ip rate constants. In order to better understand 
how the new cross-reaction feature leads to this 
dramatic improvement, we analyze which elemen¬ 
tary process immediately precedes the formation 
of an appropriate divacancy and therewith trig¬ 
gers the creation of a Pd(lOO) patch in the sim¬ 
ulations. Figure [10] shows these relative contri¬ 
butions and reveals that at lower temperatures re¬ 
duction is predominately caused by CO oxidation 
processes within the surface oxide domain, while 
at the higher temperatures the ongoing reduction 
is primarily driven by cross-reactions over y/5- 
oxide/Pd(100) domain boundaries. The latter is a 
natural consequence of the enhanced CO binding 
at the Pd(lOO) nuclei as compared to the \/5-oxide. 
At the higher temperatures CO is then predomi¬ 
nantly stabilized on the Pd(lOO) patches at the sur¬ 
face, concomitantly enhancing the cross-reaction 
compared to the regular CO oxidation within the 
surface oxide domain. This change in the mech¬ 


anism behind the oxide reduction leads to a tem¬ 
perature dependence of the reduction time that is 
impossible to grasp within a single-lattice model, 
regardless of how one tries to tweak the underlying 
rate constants. 



340 360 
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Figure 10: Relative contribution of elemen¬ 
tary processes leading to the creation of divacan¬ 
cies and therewith to the ongoing oxide reduction 
through the formation of Pd(lOO) patches at the 
surface. "Oxide reaction" denotes divacancy cre¬ 
ation through a CO oxidation reaction within the 
\/5-oxide domain, "cross-reaction" denotes diva¬ 
cancy creation through a cross-reaction over the 
\/5-oxide/Pd(100) domain boundary, and "O dif¬ 
fusion" denotes divacancy creation through diffu¬ 
sion of O atoms on the \/5-oxide. The data is ob¬ 
tained using a (20 x 20 ) simulation cell. 


4 Summary and Conclusions 

We presented a general multi-lattice Ip-kMC ap¬ 
proach that allows to describe surface morpho¬ 
logical transitions within lattice kMC, as long as 
these transitions occur between structures exhibit¬ 
ing commensurate lattices. This approach allows 
to evaluate long-time kinetic behavior subject to 
corresponding transitions and is efficient enough 
to be based on Ip rate constants. Additionally 
required compared to traditional single-lattice Ip- 
kMC is a detailed atomistic pathway for the tran¬ 
sition, which then needs to be mapped into the Ip- 
kMC model in form of new elementary processes. 

We illustrated the approach by addressing recent 
oxide reduction experiments at a Pd(lOO) model 
catalyst by Fernandes et al— Single-lattice Ip- 
kMC simulations monitoring the O coverage as 
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an oxide stability indicator can not account for 
the temperature dependence of the reduction time 
observed in the experiments. This holds, even 
when abandoning the first-principles character of 
the simulations and optimizing the rate constants 
of relevant elementary processes to fit the experi¬ 
mental data. In particular at more elevated temper¬ 
atures, the limitations in stabilizing CO at the sur¬ 
face oxide lead to reduction times that are much 
too long. This leads us to conclude that the contin¬ 
ued oxide reduction subsequent to a long induction 
period seen in the measurements is not governed 
by processes at intact oxide patches. 

Extensive DFT calculations were correspond¬ 
ingly used to establish an atomistic reduction path¬ 
way that generates local Pd(lOO) patches. An im¬ 
portant rational in this development was the real¬ 
ization that several oxide metal atoms do not nec¬ 
essarily have to change their geometric position 
during the transition from oxide to pristine metal. 
The choice to only consider pathways where cor¬ 
responding position changes are not taking place 
then leads to a drastic reduction in the space of 
possibly stable intermediates and could also be a 
useful strategy to conceive complex reconstruction 
mechanisms in other systems. 

Very good agreement with the measurements 
is obtained with the multi-lattice Ip-kMC sim¬ 
ulations on the basis of the proposed oxide re¬ 
duction pathway. The crucial new feature lead¬ 
ing to this improvement is the possibility for CO 
oxidation reactions across \/5-oxide/Pd(100) do¬ 
main boundaries. These cross-reactions constitute 
a largely varying contribution to the overall ox¬ 
ide reduction within the temperature range cov¬ 
ered by the experiments, and therewith generate a 
temperature-dependence of the reduction rate that 
cannot be effectively described by a single-lattice 
model. 

In the catalysis context an interface controlled 
reactivity is often emphasized. As long as phases 
with commensurate lattices are involved, multi¬ 
lattice Ip-kMC offers the exciting perspective to 
scrutinize this perception within comprehensive 
microkinetic modeling. The present work already 
contributes to this by highlighting that even though 
reaction barriers may be quite similar on the metal, 
the surface oxide and across the metal/oxide do¬ 
main boundary, the latter could still be the dom¬ 


inant reaction mechanism due to the steady-state 
coverages around the metal/oxide interface. For 
the present model to directly address steady-state 
CO oxidation activity primarily additional elemen¬ 
tary processes are required that describe surface 
oxidation, in particular dissociative oxygen ad¬ 
sorption at highly 0-precovered surfaces e.g. fea¬ 
turing Pd(lOO) patches surrounded by y^-oxide. 
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